Transport properties of a molecular quantum dot coupled to one-dimensional 

correlated electrons 



O 
0^ 






cd 



I 

O 

o 



> 

O 

vn 

o 
o 



X 



S. Maier^ and A. Komnik^ 

' Institut fiir Theoretische Physik, Ruprecht-Karls-Universitat Heidelberg, 

Philosophenweg 19, D-69120 Heidelberg, Germany 

(Dated: February 10, 2011) 

We analyze the transport properties of a quantum dot with a harmonic degree of freedom (Holstein 
phonon) coupled to interacting one-dimensional metallic leads. Using Tomonaga-Luttinger model 
to describe the interacting leads we construct the generating function of the full counting statistics 
(FCS) for a specific constellation of system parameters and give explicit expression for the cumulant 
generating function. In the resonant case we find the lowest order correction to the current to be 
negative and divergent when source-drain voltage approaches the phonon frequency. Via a diagram 
resummation procedure we show, that this divergencies can be repealed. On the contrary, in the 
off-resonant case the lowest order correction remains finite. This effect can be traced back to the 
strongly non-monotonic behaviour of the bare transmission coefficient (without phonon) with respect 
to the dot level energy. We calculate corrections to the noise power as well and discuss possible 
experimental implications of this phenomenon. 

PACS numbers: 73.63.Kv, 72.10.Pm, 73.23.-b 



Quantum impurities are often used to model ultrasmall 
quantum dots and contacted molecules. One of the most 
fundamental setups is the Anderson impurity model tak- 
ing into account the local Coulomb interaction^. Espe- 
cially in the case of contacted molecules, however, an 
explicit consideration of coupling to vibrational degrees 
of freedom is also desirable. This is accomplished by 
the Anderson-Holstein model^i^. Even under equilibrium 
conditions and in the absence of Coulomb repulsion it 
turns out to possess interesting properties. Recently a 
number of successful electron transport experiments car- 
ried out on contacted molecules revealed some very in- 
teresting details^"—. Probably the most prominent of 
them is the different behavior of the molecule conduc- 
tance which can grow or decline as soon as the applied 
voltage gets larger than the phonon frequencyi^i"— This 
phenomenon can be understood as follows. At zero tem- 
perature and voltages the vibrational degrees of freedom 
can be safely assumed to be frozen out and one effectively 
deals with the (noninteracting) resonant level with some 
energy A. The spectral function of the quantum dot is 
a single Lorentzian with some width T (which is related 
to the hybridization of the dot level with the electrode) 
centered around A. For the large initial transmittance 
of the system A should lie in between the chemical po- 
tentials of the contacting electrodes. On the opposite, 
for small transmittance A is below/above the chemical 
potentials. The system is virtually insulating at |A| ^ F 
because then the spectral weight around the chemical po- 
tentials position, which is necessary for transmission, is 
very small. When the phonon gets excited its spectral 
function is known to develop equidistant sidebandsi^ii^. 
The central peak at A persists but due to spectral weight 
redistribution its height diminishes. Therefore the ini- 
tially large transmission drops as soon as the vibrational 
degrees of freedom can be excited. On the contrary, 
due to the finite spectral weight in the sidebands the 



conductance grows for the out-of- resonance A. It turns 
out that in general the crossover from enhanced to sup- 
pressed transmission does not correspond to any univer- 
sal parameter constellation apart of the limiting cases of 
large/small Z^iiiii^. Nonetheless it has been observed in 
several experiments^"—. 

Thus far only noninteracting electrodes were consid- 
ered. Given the small dimensions of the correspond- 
ing devices it is very likely that the electrodes might 
in fact possess genuine one-dimensional geometry as far 
as the electronic degrees of freedom are concerned. Al- 
ternatively one might conceive a device contacted by 
e. g. armchair carbon nanotubes, which are known to 
host one-dimensional electrons. In these situations in- 
stead of conventional Fermi liquids one deals with the 
Tomonaga-Luttinger liquids (TLL). Their most promi- 
nent feature is the power-law singularity of the local den- 
sity of states in vicinity of the Fermi edge. Among other 
things it results in complete suppression of transmission 
in presence of impurities in the low energy sectors lead- 
ing to the zero bias anomaly (ZBA) ^°i^^ . As a result the 
transmission through a featureless quantum dot coupled 
to two TLLs vanishes towards small voltages and low 
temperature o^"'^^" — . The only exception is the perfect 
resonant setup when A = and hybridizations with both 
electrodes are equal to each other. Thus, contrary to the 
noninteracting electrodes, when the dot transmission can 
smoothly interpolate between perfect and zero transmis- 
sion, in the TLL setup only two low-energy transmis- 
sion regimes are possible: either zero or unity. Applying 
the above line of reasoning one might conclude that in 
the former case the current through the system starts 
to flow only after voltage gets larger than the phonon 
frequency. In the opposite case one would expect that 
the conductance of an initially perfectly transmitting dot 
would rapidly decrease beyond the threshold set by the 
phonon frequency. The goal of our paper is to under- 



stand the details of transport properties of such a setup 
and to quantify this heuristic picture. 

In order to proceed one needs a model which can 
equally good describe the off-resonant as well as perfectly 
transmitting case. While in the noninteracting case the 
resonant level model without phonon is trivially solvable 
this is not the case any more for the TLL with a generic 
interaction strength. Nonetheless, at one particular pa- 
rameter constellation the problem is exactly solvable. We 
take advantage of this and analyze the transmission prop- 
erties of the system with the phonon using the perturba- 
tion theory in the electron-phonon coupling. We model 
the system by the following Hamiltonian, 



a = (dt _|_ d)l\/2 and h = -i{S - d)/\/2 are the re- 
spective Majorana operators for the dot fermion. In the 
absence of the phonon the transport properties of the 
above system are known even on the level of the full 
counting statistics (FCS)^"— . It is accomplished by in- 
troduction of the counting field A into the Hamiltonian 
making it explicitly time-dependent due opposite sings 
of the counting field on the different Keldysh branches. 
In the model this is expressed in the replacement of 
the tunnel Hamiltonian Ht — — 1-^/276^ by the opera- 
tor T\ = —i^/2^b [^cos(A/2) — ?7sin(A/2)]. In absence of 
phonon-electron interaction the electric current is given 
by^o 



H = H\ 
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(1) 



-ffieads describes two (R/L, right/left) metallic leads in 
the Luttinger liquid state kept at different chemical po- 
tentials /ii — /Xfl = V with V being the bias voltage ap- 
plied across the dot (we use units in which e = h = ks = 
1 throughout), see for example Eqs. (2-4) of [30]. Since 
we are not interested in any spin or finite field related 
effects our dot Hamiltonian is given by 



Hdot = AdU + Qc^c + gdU{c^ + c) 



(2) 



which describes a single fermionic level with energy A 
coupled with the amplitude g to the local oscillator with 
frequency O. d, d^ and c, c^ are the respective annihi- 
lation/creation operators for the electron and oscillator, 
^tunn is responsible for the particle exchange between 
the electrodes and the dot, 
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where jr^l are the (constant) tunneling amplitudes. The 
tunneling is assumed to be taking place locally at x = 
in the coordinates of the respective electrode. The last 
term is the capacitive coupling of the electrodes to the 
dot 
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The present problem allows for an exact solution for g = 
and Tomonaga-Luttinger interaction parameter K ~ 
1/2 while Uc — 2ttvf- The emergent structure is that 
of the Majorana resonant level model (MRLM) extended 
by the vibrational degree of freedom^^, in the case of the 
symmetric coupling 7l — 7i? = 7 

H = Ho[^, T]]-i[A + g{c'' + c)] ab + t-fb^{Q) + nc^fc . (5) 

The first part of the transformed Hamiltonian describes 
the electrodes in terms of new Majorana fermions ^ and 
V, 



Ho = i dx 



il{x)d^il{x) + i{x)d,i{x) + Vi{x)T^{x) .(6) 



I = Go du {fiL - nR)D{uj) , 



(7) 



where Gq = 2e^//i is the conductance quantum, riL^n are 
the Fermi distribution functions in the respective elec- 
trode and 
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is an effective transmission coefficient for the new 
fermions. Its nonmonotonic behavior around a; = for 
zero and finite A is the reason for the full transmission 
suppression in the off-resonant case22r24. The system un- 
der consideration is not the only one with the effective 
transmission coefficient of the form ^. It is also encoun- 
tered in the exact analytic solution of the two-terminal 
Kondo model at the Toulouse point^i as well as in the 
transmission coefficient for the physical electrons in the 
parallel double quantum dot setup^>^. The generaliza- 
tion of the phenomena discussed below to these setups is 
thus straightforward. 

The cumulant generating function (CGF) for the sys- 
tem at g ^ factorizes as \(X) = TCo(\)Xn(^), where 
XoW is the CGF calculated in |27H29l |. The correction 
due to electron-phonon coupling is found from 



InXg(A) = L-9lc'ir[cHr)+c{T)]a{T)b(r)\ 



(9) 



where the expectation value is taken with respect to the 
operator Ho [£,, rj] — z A a5 -I- flc^c + T\ . The time integra- 
tion runs along the full Keldysh contour C and the super- 
scripts k,l — zL refer to the backward/forward propagat- 
ing branches. In the leading order perturbation theory 
there are two different contributions: (i) the one gener- 
ated by the 'tadpole' diagram 
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(ii) the one given by the 'shell' diagram 
^ ,, ,__■_ Jc 



k,l = ± 



Dl (ri2) D[% (-ri2) - i?,^i (n^) D^ {-r. 



(11) 



Here, we introduced the phonon Keldysh Green's 
function (GF) Gph{Ti2) = -« (ic c(ti) c^ (rs)), the 
Keldysh functions of the quantum dot Dap {T12) — 
— i (Tc a (ti) /3 {t2));^ where a and f3 take values in the 
set {a, 6} and T12 = ti — T2. At this point it is convenient 
to transform to energy variables 
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and introduce generalized (Keldysh) polarization loops 
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In the resonant case (A — 0) the mixed dot Keldysh GFs 
(a 7^ /3) vanish identically and Daa (i-^) becomes diagonal. 
Hence, we are left with the calculation of the (anti-)time- 
ordered components of the polarization loop. In the zero 
temperature limit the first order correction to the full 
counting statistics can be calculated exactly. One finds 



InxW 



n- 2 e* 2 r tan ^ ( ^e ' 



"tanh-i(^) 
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Thus the correction to the transport current has the fol- 
lowing form, 



j(i) 
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tan-i (V/T) - 2rO tanh"^ {V/fl) 
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Expressions for noise and higher cumulants can easily be 
determined by taking higher order derivatives of the CGF 
with respect to A and setting it to zero afterwards. The 
calculation of the correction to the full statistics at finite 
temperature is quite involved. It is, in fact, more con- 
venient to set up explicit expressions for the individual 
cumulants before the energy integration. A remarkable 
property is the absence of any sharp energy threshold 
usually found in phonon affected transport between un- 
correlated electrodes. Those are usually attributed to 
the onset of inelastic processeaii"— , when individual elec- 
trons may loose/gain energy n during tunneling between 
the leads a.tV>fl. In the present case of MRLM the sit- 
uation is different. The Majorana fermions describe the 
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FIG. 1: Main graph: first order current correction at finite 
temperature. The parameters are n/F = 1 and T/T = 
0.1,0.05,0.01 (solid, dashed , dashed-dotted) . Inset: first 
order current corrections in case of zero temperature for 
r/fl = 1,2,3 (solid, dashed, dashed-dotted) 



collective excitations of the TLL state in the leads (kinks 
and antikinks), rather then individual electrons. That 
is why even at T = there is no sharp energy thresh- 
old. Another feature is the logarithmic divergency for 
F — > r2, which does not survive at finite temperatures, 
see Fig.[T] It can be attributed to the non-monotonic be- 
havior of the effective transmission coefficient (|8]) at zero 
and finite off-set A. The transmission for the quasipar- 
ticles drops immediately to zero as soon as A becomes 
off- resonant. This is exactly what occurs when the har- 
monic degree of freedom is excited - the coupling g gen- 
erates an effective non-zero offset. In this voltage regime 
the problem becomes non-perturbative in the electron- 
phonon coupling and therefore requires analysis by more 
advanced methods. One possible route would be an RPA- 
like diagram resummation. To perform this program, we 
express the CGF in terms of the adiabatic potential^!, 
which itself can be rewritten as a function of the exact 
Dhb Keldysh GF only. It is approximated by the RPA- 
type summation (see Fig. [5] for the diagrammatic struc- 
ture). Interestingly, within this approximation one finds 
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FIG. 2: Summation of an infinite series of diagramms. Dtt, 
Daa represents the Keldyshfunction of the b resp. a Majorana 
fermions as defined before with A set so zero and the wiggly 
lines represent the phonon propagators 



for the CGF the Levitov-Lesovik formulai^ 

lnx«P^ (A) = r I g In [1 + Df (c.) (e^^ - l) {n, - nn) 

(15) 
with the effective transmission coefScient 



Dfiuj) 



r2 (^2 _ ^2^ 



g4a;2 - 252^2 (^2 _ ^2) j^ (^2 + r2) (^2 _ Q2y 
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which has a three maxima (perfect transmission) at w = 
0, zLy^g'^ + $72 and minima (perfect reflection) at a; = 
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FIG. 3: Main graph: Full current in RPA approximation. The 
parameters are Q./T = 1 and g/T = 0.25, 0.5, 0.75, 1 (solid, 
dashed , dashed-dotted, dotted lines, respectively). Inset; 
the current correction in RPA approximation Q,/V = 1 and 
g/T = 0.25, 0.5, 0.75, 1 (solid, dashed, dashed-dotted and dot- 
ted lines) 

In Figs. [3] and [5] the current and noise in our approx- 
imation are depicted. One observes that for small volt- 
ages the current increases nearly linear. For voltages near 
n one finds a plateau, e. g. the current enhancement 
is suppressed by the electron-phonon interaction. This 
feature does not occur exactly a,t V = fl (see Fig. [3]). 
One finds that the maximal reduction of the current is at 

V ^n^ 



1 



2fP 



or if one assumes Q^ g, V ~ n + jj^. 

This kind of a shift by ^ is one normally produced by a 
polaron (Lang-Firsov) transformation^^. 

The behavior of the shot noise turns out to be even 
more interesting, see Fig. |4l In case of zero temperature, 
to the lowest order in g one again finds a log-divergent 
contribution for voltages approaching the phonon fre- 
quency. However, for a special configuration F = fi the 
singularity cancels out and the shot noise interpolates in 
a regular way between small and large voltage limits. For 
O > F the correction to the shot noise develops a new 
feature: it changes sign, see Fig. 2] inset. Just as in the 
case of the transport current, the resummation heals the 
singularity, see Fig. [5] 
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FIG. 4: Main graph: g^ correction to the shot noise at fi- 
nite temperatures. The plotting parameters are Q/T = 1 and 
T/r = 0.1,0.05,0.01,0 (solid, dashed, dashed-dotted, dot- 
ted lines respectively). Inset: g^ noise corrections in case of 
zero temperature. The plotting parameters are F/fi = 1,2,3 
(solid, dashed, dashed-dotted lines) 
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FIG. 5: Main graph: Shot noise in RPA approximation. 
The parameters are Q./V = 1 and g/T = 0.25, 0.5, 0.75, 1 
(solid, dashed , dashed-dotted and dotted lines). Inset: 
noise correction in RPA approximation for Q,/V — 1 and 
g/r — 0.25, 0.5, 0.75, 1 (solid, dashed, dashed-dotted and dot- 
ted lines). 



There is an important difference in the behavior of the 
noise and current corrections in the high voltage limit. 
While the latter tends to vanish the former remains finite. 
It can be understood as follows. As already discussed in 
the introduction, for the current only the overall spec- 
tral density between the chemical potentials in the leads 
is essential. For V ^ il the phonon excitations are less 
effective in squeezing the spectral weight beyond the volt- 
age window, so the current should approach the unitary 
value already in the zeroth order in g. For the noise, how- 
ever, the rate of phonon (de)excitation is important. In 
the limit y — > oo it achieves its maximal value, thereby 
generating a finite contribution seen in Fig. 3] However, 
this feature isn't observed after the diagram subsumma- 



tion procedure. 

We now turn to the ofF-resonant system A 7^ 0. It is 
again convenient to perform the A differentiations first 
and integrate over energy afterwards. Here two different 
regimes of weak |A| < r/2 and strong detuning |A| > 
r/2 emerge,22*26 

In the foUowing we consider the case of weak detuning. 
Then, the tadpole contribution yields 
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where we have defined f2± = r/2 ± ^/ {V^ / A - A^) . The 
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FIG. 6: Tadpole correction to the current. The parameters of 
the plots are Q./r = 10 and A/F takes the values 0.25, 0.2, 0.17 
(solid, dashed and dot-dashed lines, respectively) 

corresponding current correction is plotted in Fig.|6l Just 
as in the case of transport between two uncorrelated 
electrodesi^, its simple functional dependence on the 
phonon frequency hints at its dominance in the (Born- 
Oppenheimer) limit of the very 'slow' phonon^^i^. The 
parameter F mainly determines the width of the peak 
and A its position. Because of the simple form of this 
contribution, we will omit it in further illustrations. 

The contributions from the shell diagram are more so- 
phisticated: 



/ii:ln = ^/_'|bi(y)+. 



J2{y)+j3{y)], (18) 



where we have 
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Do [y ± n) 



A^){y±n) {y + ny^A^^ 
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y{y^ 



Dl{y) 



where Do (y) — (y^ — A^) + y^T^. The expression von 
J3 is quite lengthy, see Appendix for details. Similar to 
the tadpole contribution this part is also always posi- 
tive. In contrast to the resonant case, the corrections re- 
main hnite for every voltage. For V close to the phonon 
frequency we observe a double step like feature instead. 
This can be attributed to the double-maximum behav- 
ior of the transmission coefficient ([5]), the denominator 
Do{y) of which enters all relevant Green's functions and 
thereby affects the voltage dependence. The linear scal- 
ing of the distance between the dips is shown in Fig. [T] 
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FIG. 7: Main graph: shell diagram correction to the cur- 
rent. The parameters of the plots are A/F — 1 and Q/T 
takes the values 9, 10, 11 (solid, dashed, dashed-dotted). In- 
set: double step like feature for the parameters Q/T — 10 and 
A/F = 1,1.5,2 (solid, dashed, dashed-dotted). In both, the 
main graph and the inset, we chose strong detuning, because 
in this limit electron-phonon interaction effects are more pro- 
nounced. 

To conclude, we investigated the interacting resonant 
level model in presence of a harmonic degree of freedom 
coupled to the quantum dot. We observe that in the res- 
onant case, where the system is initially perfectly trans- 
mitting, finite electron-phonon coupling leads to negative 
corrections to the current. In the zero temperature limit 
we identified a strongly non-perturbative regime where 
the current correction is log-divergent and performed an 
RPA-like resummation of divergent diagram contribu- 
tions, which turned out to produce a plateau- like feature 
in the full current- voltage characteristics of the system. 
We believe that this behaviour is generic in all setups 
with TLL electrodes also beyond the chosen parameter 
constellation. Single wall carbon nanotubes (SWNTs) 
are known to be typical realizations of the TLL electronic 
statsi^"— Therefore we expect the above strong conduc- 
tance suppression phenomenon to be observable in exper- 



iments on molecular quantum dots coupled to SWNTs. 
In the opposite off-resonant case, when the system with- 
out the phonon has zero conductance, we observe conduc- 
tance enhancement due to electron-phonon interaction. 
For voltages comparable to phonon frequency wc find a 
double-step like feature in the lowest order perturbation 
expansion in electron-phonon coupling. Contrary to the 
resonant case no singularities are observed. 
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Appendix A: Tadpole contribution in case of strong 
detuning 



In case of strong detuning A^ > ^ one finds for the 
static contribution the expression 
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Appendix B: Explicit expression of the ja current 
contribution 



With the definitions 
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one finds for J3 
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